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Abstract. CsCuCla is a ferromagnetically stacked triangular spin-1/2 antiferromagnet. We discuss models 
for its zero-temperature magnetization process. The models range from three antiferromagnetically coupled 
ferromagnetic chains to the full three-dimensional situation. The situation with spin-1/2 is treated by 
expansions around the Ising limit and exact diagonalization. Further, weak-coupling perturbation theory 
is used mainly for three coupled chains which are also investigated numerically using the density-matrix 
renormalization group technique. We find that already the three-chain model gives rise to the plateau-like 
feature at one third of the saturation magnetization which is observed in magnetization experiments on 
CsCuCb for a magnetic field perpendicular to the crystal axis. For a magnetic field parallel to the crystal 
axis, a jump is observed in the experimental magnetization curve in the region of again about one third of 
the saturation magnetization. In contrast to earlier spinwave computations, we do not find any evidence 
for such a jump with the model in the appropriate parameter region. 



PACS. 75.10.Jm Quantized spin models - 
- 75.50.Ee Antiferromagnetics 



75.45.-l-j Macroscopic quantum phenomena in magnetic systems 



1 Introduction 



CsCuCla is in several respects a quite unusual stacked 
triangular (anti)ferroniagnet (see e.g. [y for a review of 
the subject). Among others, three-dimensional features 
are vital in this compound: The ferromagnetic stacking is 
just an order of magnitude larger than the antiferromag- 
netic in-plane coupling while in most other materials the 
coupling constants differ at least by two orders of magni- 
tude. Furthermore, the smaller antiferromagnetic coupling 
is crucial since it gives rise to a non-trivial behaviour of 
CsCuCla in the presence of an external magnetic field. 

The behaviour of CsCuCls in strong external magnetic 
fields and at low temperatures has been studied already 
some time ago (see e.g. g]). One observes different be- 
haviour, depending on whether the magnetic field is ap- 
plied along or perpendicular to the c (crystal)-axis. For a 
magnetic field along the c-axis, there is a jump in the mag- 
netization curve at values of the magnetization of about 
one third of the saturation value. In contrast, one finds a 
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plateau in the same area if the magnetic field is applied 
perpendicular to the c-axis|^. 

The spin in CsCuCla is carried by Cu^+ ions and thus 
is a spin 1/2. So, a proper treatment of quantum fluc- 
tuations is important. However, in order to be able to 
treat the three-dimensional situation, most theoretical ap- 
proaches have so far been either phenomenological P,0Jq| 
or used first-order spinwave theory |0,g,||,|o[|^ . 

The goal of the present paper is to treat directly a 
spin-1/2 model with realistic XXZ-anisotropies. We will 
use series expansions around the Ising limit, numerical di- 
agonalization and degenerate second-order weak-coupling 
perturbation theory. 

When one solves a classical or Ising-version of this 
model, it becomes clear that essential features are already 
captured by a three-chain variant. The bulk of our paper 
will therefore concentrate on three antiferromagnetically 
coupled ferromagnetic chains. The simplification to a one- 
dimensional situation also opens the possibility of apply- 
ing methods which work particularly well in one dimension 
like the density-matrix renormalization group (DMRG) 
method ||ll|,0. 



^ Actually, a plateau with (M) = 1/3 can also be observed in 
other triangular antiferromagnets (see e.g. [H] for recent rather 
clear examples). 



A. Honecker et al.: A Spin-^ Model for CsCuCls in an External Magnetic Field 



The purpose of the present paper is to gain some qual- 
itative insight and demonstrate that new theoretical ap- 
proaches to CsCuCls are viable. It should be possible to 
develop all these approaches further to improve the accu- 
racy of the results, to obtain a more realistic modeling or 
to compute other quantities. 



2 The model 

Now we proceed to give a more accurate description of our 
model and the parameter region we are interested in. We 
model CsCuCls by a stack ed t riangular lattice antiferro- 
magnet as sketched in Fig. 2.1. 




Fig. 2.1. The stacked triangular (anti)ferroniagnet in a mag- 
netic field h. The bold lines denote our three-chain model. 

The quantum-spin Hamiltonian corresponding to Fig. 



2.1 is given by (compare also |8| 



N L 



^ ~ '^ A^A^ ^'^^i,xSi,x+l + 2 \^i,x^i,x+l + ^'i,x^i,x+l) 



= 1 x = l 



'^J 2^ 2^ ^'2^i,xS'],x + 2 \^i,x^j,x + ^i,x^],x) 

-hY^Sl,, (2.1) 

i,x 

where the Si^x are spin-1/2 operators and ft, is a dimen- 
sionless magnetic field. The notation (j,j) denotes neigh- 
bouring pairs of the N chains. Here, the number of chains 
N will be kept fixed while we are interested in the limit 
of infinite L. Mostly, we choose Ai = A2 and then denote 
both of them just by 'Z\'. 

Quite often, CsCuCls is modeled by a Hamiltonian 
which includes a Dzhalozhinski-Moriya interaction with a 
vector pointing along the ferromagnetic chains in order to 
account for a modulation in the magnetic structure with 
a period of slightly more than 70 layers along the c-axis 
p3| . However, by a local unitary transformation this can 
be traded for a contribution to the anisotropy A2 and a 
surface term which is then discarded (see e.g. ||]). Further 
contributions to the anisotropics come from crystal fields. 



In this way, the z-direction in spin space is identified with 
the direction along the ferromagnetic chains (the c-axis) 
in real space. 

J will be ferromagnetic (J < 0) and J' antiferromag- 
netic (J' > 0). In order to describe the material CsCuCla 
one uses J' w -J/6 and A slightly less than one (l|,|l|,|l6|,0 
The actual values used should not matter so much; it 
should however be noted that CsCuCls is described by 
weakly coupled ferromagnetic chains ~ though not so weakly 
coupled as is the case for many other triangular lattice an- 
tiferromagnets 0]. A Monte-Carlo simulation of an XY~ 
model takes J' = —J/ 10 [E8| showing the influence the 
weaker coupling has on the critical behaviour of triangu- 
lar antiferromagnets. To make the situation clearer in the 
study of the few-chain problem we choose a bit stronger 
interchain coupling, e.g. J' = —J/3 to compensate par- 
tially for the missing neighbouring chains. 

The (i, j)-summation in (2.1) sh ould run over a tri- 
angular lattice as sketched in Fig. 2.1. However, since 
J' <C I J|, a reasonable first approximation should be ob- 
tained by retaining only one site as a representative for 
each of the three sublattices of the triangular lattice. This 
leads to the model of three coupled chains which is indi- 



cated by the bold lines in Fig. 2.1 



Since we are interested in the behaviour of CsCuCls at 
low temperatures, we actually simplify to zero tempera- 
ture in the present paper. We further restrict to magnetic 
fields applied along the anisotropy axis. The magnetiza- 
tion 



IM) = 



N L 



(2.2) 



\i=l x=l I 

is then given as the expectation value of a conserved op- 
erator. This is technically useful since it permits one to 
relate all quantities in a magnetic field h to those with 
a given magnetization (M) at zero field. In the presence 
of an anisotropy Z\ 7^ 1, the effect of a magnetic field 
along the z-axis is different from a magnetic field in the 
xy-plane. Because of the aforementioned simplification we 
consider only the former case in the present paper. 

Possible plateaux in the magnetization curve can be 
most easily read off in the strong-coupling limit J' ^ | Jl. 
For a fixed number N of coupled chains one finds [[l9| , P0| 
that for J = the only possible values of the magne- 
tization are (M) = -1, -1 -f 2/iV, ..., 1 - 2/iV, 1. In 
particular, for all N that are divisible by three, a plateau 
with (M) = 1/3 is possible (the simplest case is iV = 3). 
However, this strong-coupling argument does not ensure 
that a certain plateau actually does occur for given val- 
ues of the parameters - this issue requires a computation 
of the magnetization curve (at least in the vicinity of the 
plateau- value of the magnetization) . 

Before we proceed with the discussion of the magne- 
tization process, we note that a straightforward compu- 
tation yields the upper critical field for A^ = 3 cyclically 
coupled chains 



Kc^J'{A2 + ^) + J{A^-l) (2.3) 

for J' > and J < 0. This result is based on the as- 
sumption that the transition proceeds with a single spin 
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flip and is therefore valid for Z\i small enough (in partic- 
ular A-\ < 1). If the transition to saturation becomes first 
order, ( ^.3[ ) ceases to be valid. 



3 Ising expansions for three chains 



The Ising limit of iV = 3 ferromagnetic chains (2.1) which 



are mutually antiferromagnetically coupled is easily un- 
derstood: At zero temperature, only states with {M) = 
±1 and (M) — ±1/3 are realized. The groundstate at 
(M ) = 1/3 is given by a configuration in which all spins 
in two chains point up and those in the third chain point 
down. Perturbation expansions in Z\~^ around this state 
and the ones with a single spin fiipped respectively to it 
are readily performed. One finds the lowest gaps for the 
single spin excitations with states that are translational 
invariant along the chain direction, and antisymmetrized 
between the two equivalent chains for an additional spin 
flipped down. 

With this information one proceeds in the same way 
as for the triangular lattice |gl[ (the precise method has 
been summarized in |£2| ) to flnd sixth-order series for the 
gap of a single flipped spin. The series for general J and 
J' can be found in appendix A. Here we just present the 
specialization to J' = —J/3: 



\J\ +6 



11 



1993 



5843545 



"'C2 
1^1 



432Z\2 
202825209149 
"^ 2824576634880/4'' 
4A 1 11 

~3" ^"^^ mZ ^ 432/42 



1558903 



18289152Z\4 224682232320/45 



72576Z\3 109734912Z\4 

+ 0{A-'^), (3.1) 

(3.2) 
0{A-^). 



3887 
54432Z\3 
33098775077 



It should be noted that for sufficiently large A it is actually 
favourable to flip large clusters of spins rather than a single 
one if one changes the applied magnetic fl eld. There fore , 
the tran sition at la rge A is flrst order and (3.1) and (3.2) 
(or (A.l) and ( |A.2| )) are not the boundaries of the (A/) — 
1/3 plateau. Nevertheless, in contrast to other Ising series 
such as those for the triangular lattice [^ , these series are 
remarkably well behaved down into the region Z\ « 1 for 
the values of parameters we are interested in (J' small). 
The ending point of the (M) = 1/3 plateau is determined 
by the condition hc2 (Ac) — ft-ci (^c) = 0. Fo r exa mpl e, for 



J'/l J| = 1/3 we can use the raw series (3.1) and (3.2) and 
flnd that the plateau ends at A^ = 0.965, 0.876 or 0.997 
if one truncates the series at fourth, fifth or sixth order, 



respectivel y. At J' /\J 



= 1/6 convergence is a bit better 



and using (A.l) and (A. 2) one finds the ending point at 
Z\c = 0.972, 0.954 and 0.985 for fourth, fifth and sixth 
order, respectively. 

This already indicates that the precise value of the 
XXZ-anisotropy is crucial for explaining the presence or 
absence of a plateau with (M) = 1/3. 



4 Exact diagonalization for three chains 

Next we present exact diagonalization results for three 
cylindrically coupled chains with length L = 8. Fig. 4.1 
shows the magnetic phase diagram for the su(2)-synimetric 
situation A — 1. The lines show the magnetic fields at 
which a transition occurs between the states with the 
magnetization indicated in the figure (for L = 8, only 
(M) — to/12 with 771 = —12, . . . , 12 can be realized). One 
observes a plateau with (M) — 1/3 for all J' > 0. The 
transition to saturation clearl y is given by h^c = §J', 
as it should be according to (|2.3| ). In the ferromagnetic 
regime J' < 0, the magnetization jumps from (M) = — 1 
to (A/) = 1 as ft- passes through zero. 

We would like to remark the striking smoothness of 
the transition lines in the weak-coupling regime J' <^ \J\ 
in Fig. 4.1. This suggests that expansions in J' should be 
a viable way to understand the magnetization process of 
three coupled chains (at least in a situation with some 
su{2) symmetry preserved). This will be the subject of a 
later section. 

For J' 3> \J\, we can use the effective Hamiltonian 
given in |19| to determine the lower boundary of the {M) — 
1/3 plateau. By exact diagonalization we find with A = 1 
that ftci/l J| = 0.42096, 0.41880, 0.41779, 0.41723, 0.41689 
for L = 8, 10, 12, 14, 16, respectively, and in the limit J/ J' - 
0. By comparison with Fig. ^J we see that the lower 
boundary of the (Af) = 1/3 plateau must therefore still 
increase beyond the right boundary of the figure. 



1.25 



h/IJI 0.75 




0.25 - 



Fig. 4.1. Magnetic phase diagram for three coupled chains 
with L = 8 and A = 1. 



Now we fix the coupling constants to J' = | Jl/3 and 
look at the dependence on the XXZ-anisotropy A. The 
resu lting magnetic phase diagram for L = 8 is shown in 
Fig. IJ. 

One observes that for /4 > 1.14 only the values (AI) = 
±1/3, ±1 occur. The corresponding plateaux are sepa- 
rated by first-order transitions in this region. Outside this 
region, the lower boundary of the (Af ) = 1/3 plateau and 
the transition to saturation seem to remain second-order 
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0.4 --_5M> = 2/3 

h/lJl 0.3 ,__ <M>^r77J 
<M> 



Fig. 4.2. Magnetic phase diagram for three coupled chains 
with L = 8 and J' — JI/3 . The dashed lines denote the sixth- 
order series ( |3.l| ) and (3.2). The diamonds denote extrapolated 
values of the boundary of the (M) = 1/3 plateau (see next 
section). 



transitions. The transition to saturation is given by ( p.3[ ) 
in the region where it is of second order. 

The (M) ~ 1/3 plateau can be observed in Fig. [4.2| 
for Z\>0.98, i.e. it survives a small amount of XY-like 
anisotropy. The transition at its upper boundary becomes 
first order for weak Ising-like anisotropics. For L = 8 
flipping two spins is favoured over flipping just one for 
Z\>1.05. A guess about what happens in the thermody- 
namic limit L — > oo is not obvious. Furthermore, in the 
vicinity of Z\ ~ 1.13 there is some evidence for a jump in 
the magnetization curve at /i = 0. 

The dashed lines in Fig. |4.2| show the series expansions 

(3.1) and (3_^) in A^^. They denote the energy-cost to 
flip a single spin above the (M) = 1/3 plateau and should 
therefore be compared to the corresponding full lines. The 
agreement of these sixth-order series with the numerical 
data is quite good in view of the fact that the window 
shown in the figure is rather far from the regime A — > 
cx) where the series were derived. These series point to a 
closing of the plateau in about the same region indicated 
by the numerical data. 

We conclude this section by noting that the transition 
to full magnetization is compatible with the Dzhaparidze- 
Nersesyan-Pokrovsky-Talapov (DN-PT) universality class 
p3| , p^ . From a technical point of view, it is important that 
conservation of the magnetization M reduces the dimen- 
sion of the space of states considerably such that we can 
apply the exact diagonalization technique to large lattices. 
Actually, for many purposes it is sufficient to consider only 
two-spin deviations p5| which can be easily treated nu- 
merically for volumes of around 1000 spins. Nevertheless, 
for J < one observes much stronger crossover effects 
than for J > |20l. 



5DMRG 

Because of the one-dimensional topology of the three-leg 
ladder we have used DMRG [ PLl] (for a detailed review of 
the DMRG method see also ||12|| ) in order to determine 
more accurately the value of the anisotropy A where the 
plateau in the magnetization occurs. Similar DMRG com- 
putations for the magnetization process of a three-chain 
model with J > have been carried out in |2^,^ (see also 
P8[ for a DMRG study of other spin ladders in a magnetic 
field) . Details on our implementation of the DMRG pro- 
cedure can be found in appendix B. 

To do the calculations, th ree spins Si^x, ? = 1,2,3 for 
fixed X (one triangle in Fig. 



2.1 



were combined to one 
site. In this way one ends up with a spin chain of L sites. 
We have calculated the (M) = 1/3 plateau for ladders of 
90 spins (30 sites in the corresponding chain model) for 
different values of the anisotropy A. 



0.35 




Fig. 5.1. Upper and lower critical fields as a function of 1/L 
for different values of A. Filled symbols denote DMRG results 
and open symbols data from exact diagonalization. 



Fig. |5.l| shows the upper and lower critical fields as 
a function of the inverse chain length 1/L. For large A 
(> 0.98) the (M) = 1/3 plateau also exists in the thermo- 
dynamic limit. In the case of small A (< 0.935) the data 
is compatible with a vanishing plateaun. For A < 0.8 the 
plateau width follows a 1/L law to a good approximation 
(this form is a strong indication of a vanishing plateau 
width). 

The interesting point is the threshold where the pla- 
teau opens. However, the infinite-system algorithm of the 
DMRG method fails already for short chain lengths when 
A approaches the region around 0.95 where the threshold 
is expected. The reason for this is not yet clear. Therefore 
in Fig. 5.1 results from exact diagonalization are plotted 



Variation of the truncation level in the DMRG procedure 
indicates that the data may have errors up to 10^'^ J . . . 10~^ J 
for the larger system sizes. 
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for these values of A. The finite-size data for the plateau- 
boundaries has been extrapolated to infinite L using ei- 
ther a polynomial fit in 1/L or a fit with corrections of 
the form exp(— cL). The former fit should work well for 
small A while the latter one is justified when a plateau 
is present, i.e. for the larger values of A. The result s of 



the extrapolation are shown by the diamonds in Fig. 4.2 



The error of the extrapolation is smaller than the size of 
the symbols, but there are additional systematic errors 
(e.g. in the DMRG truncation procedure) which may be 
somewhat larger. Therefore, all the data taken together 
indicates that the plateau closes indeed at Zi ~ 0.95 - 
at least it becomes unobservable in this region and for 
all practical purposes can be considered to be absent for 
A < 0.95. 



6 Weak-coupling expansions 

Since the coupling constants for CsCuCla lie in the weak- 
coupling regime J' ^ |J|, it is natural to try to analyze 
( |2.l[ ) by an expansion in J' around the decoupled point 
J' = 0. Suc h an approach is even more strongly suggested 
by Fig. 4.1 which shows that at least for A^ = 3 and A — 
the transition lines are very smooth functions of J' /\J 



1, 



6.1 First order 

We choose Z\i = 1 in (^^) for the weak-coupling expan- 
sions because of two reasons {A2 will however be retained 
as a parameter). Firstly, we can then exploit the su{2) 
symmetry of the decoupled chains in order to immediately 
write down their groundstates for a given magnetization. 
Secondly, in this case we can expect a smooth dependence 
on the coupling J' (see Fig. 4.4), which is not present 
otherwise (this is ind icat ed by numerical diagonalization 
- compare also Fig. 4.2). A consequence of this choice 
Al = 4 is that a given magnetization can be arbitrar- 
ily distributed among the individual chains. This means 
that the groundstate is highly degenerate and one has to 
perform degenerate perturbation theory. 

The groundstate space of a single su(2)-symmetric fer- 
romagnetic fleisenberg chain of length L is the spin-L/2 
representation. Thus, one can immediately write down the 
groundstate with a given total S" eigenvalue m: 



\j^L/2,m) 



1 



V^ 



(s- 



\jj), 



(6.4) 



where <S~ is the total step operator for a single chain. 
Here j denotes this spin and m the eigenvalue of the z- 
compone nt o f the total spin operator. The normalization 



A/j,m in (3.4) should be chosen such that all states \j,m) 



are normalized to one. Otherwise it is sufficient to know 



that the number of terms on the r.h.s. of (O) is ( _■ 



2j 



Because of translational invariance one readily finds 



that 



1 772 777- 

(j>'7i|S'f |j,m) = -{j,m\S''\j,m) = — = — . (6.2) 
L L 2j 



A simple combinatorial consideration further leads to 



{j,m+l\S:^Jj,m) 



ViJ -m){j + m+ 4) 
2j 



/■ lie- I- \ V(j-"i-h4)(j-hm) 
(j,m-4|S'^_^|j,m) = — . (6.3) 

Note that due to translational invariance the matrix ele- 



ments ( 6^ ) and ( 6^ ) do not depend on the position of the 
spin-operator x. The summation over x in ( ^.4[ ) therefore 
just yi elds a fact or L = 2j. 

In (6.2) and (6.3) one recognizes the matrix elements 
of the J operators of a spin-j sM(2)-representation nor- 
malized by a fact or 4 /(2j). Thus, the interaction part of 
the Hamiltonian (|2.4D can be rewritten in first order as 



.JkJf 



JkJi 



hY^Ji 



■' (k.l) k 

(6.4) 
This form is particularly useful for an analytical treat- 
ment. Firstly, one sees that for L ~ 2j ^ 00, the first- 
order interaction maps to a problem of classical spins 
where each spin stands for one chain. In the groundstate 
of a classical spin model on the triangular lattice, all spins 
on one sublattice point in the same direction. Therefore, a 
triangular arrangement of ferromagnetic chains and three 
coupled chains become equivalent at first order in J' - 
only the coupling constant is rescaled by the different co- 
ordination number. 

Using this mapping to a classical spin model, a number 
of conclusions can be drawn which apply to three coupled 
chains as well as the full triangular arrangement: 



The entire magnetization curve is linear at A2 = 4. 
For A2 > 4, the magnetization jumps from {Al) ~ 
—mo to (M) = niQ with mo = {A2 — 1)/ {3(A2 + 4)) as 
h passes through zero (see also Fig. [4.2[ ). 
For A2 > 1 there is a plateau with (M) — 4/3. Its 
lower boundary is given by ft-d ~ zJ' /A where z is the 
number of nearest neighbours (z = 2 for three coupled 
chains and z = 6 for an underlying triangular lattice). 
The upper boundary is more difficult to determine an- 
alytically - a lower bounc^ is given by hc^ > ZA2J' /A. 
The asymptotic behaviour of (M) as a function of h is 
linear close to saturation or to the boundaries of the 
(M) = 4/3 plateau: 



(M) - Mc (X h ~ he 



(6.5) 



With N = 3 chains, the first conclusion can actually be 
obtained for finite L using a different argument which is 
also useful for determining the groundstate degeneracy: 
For A2 ^ 1 and A^ = 3 one can further rewrite eq. (^^) 
as 



«'-J(J. 



,,^,,f-i!:ii±ii 



-h(j; + ji + jS) 



(6.6) 



The important point is that hc2 > hci for A2 > 1. 
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This shows that the energy is a quadratic function of the 
total spin (or the magnetization (M)). Therefore, all steps 
in the finite-size magnetization curve acquire equal width 
at Z\2 = 1. Also the groundstate degeneracy can be easily 
inferred from the representation (p]l). It is found to be 
maximal at (M) = 1/3 where it is equal to i + 1 and 
minimal at (M) = 0, 1 where the groundstate is unique. 

The treatment of the second order is going to be based 
on a numerical determination of the groundstate of the in- 
teraction Hi at TV = 3 and finite L. It is then important 
that the groundstate should be non-degenerate at first or- 
der. This precludes an analysis of the second order for the 
case A2 — I where after treatment of the first order one is 
still left with large degeneracies. For A2 y^ 1, the ground- 
states with a given magnetization are at most twofold de- 
generateQ. This twofold degeneracy is due to the spatial 
symmetries of the three-chain model which can be used 
to lift it completely (alternatively, such a two-fold degen- 
eracy can be ignored since it does not affect the results). 
Thus, in particular the case we are mainly interested in, 
namely A2 < 1 is amenable to a numerical treatment up 
to the second order. 



6.2 Second order 



The fc 7^ generalization of (|6.3|) is found to be given by 



{j,m + l,k\S:^Jj,m) = 



(j,m- l,k\S,Jj,m) 



^—ikx 


1 U - m - l){j - m) 


2j 


V 2j 1 


■ 


/(j+m-l)(j + m) 



2j 



2j-l 



(6.10) 



With the matrix elements (6^) and ( |6.10 ) one can now 
easily write down the matrix elements for the interaction 
between two chains. It should be noted that due to trans- 
lational invariance, spinwave states with momenta k and 
—k have to come in pairs (in other words: ^^ g-»(fe-(-fc )x _ 
L6h -^k' n)- Thi s lea ds to a cancellation of the phase factors 
in (6^) and ( 6.10 ) and thus makes the matrix elements 
of the interaction fc-independent. Then the fc-summation 
can be pulled out from the second-order matrix. So, one 
obtains for the second-order contribution 



(,',iv("" 



V 4 



Ho) F|Vi) = £~\i^i\HlHji\tl,i) , 



(6.11) 

where |V'i) is the eigenstate of the first-order matrix Hj. 
The energy denominator is given by 



If one goes beyond the first order, one has to consider also 
excitations that are created by the interaction from the 



ferromagnetic groundstates (6.1). At second order, only a 



single spin can be flipped in each chain. To deal with the 
second order perturbation it is therefore sufficient to look 
only at spinwave states 



£-'^ 



1 ^-^ 1 

2|J| ^ 1-cos 



{L-l){L + l) 



iJl^^l-cos(^) 12|J| 

(6.12) 
In order to write Hn in a compact form, we introduce 
operators Ef , Ef acting in the groundstate space of the 
ith chain as 



I j, m, k) 



1 



''^j.m ^ — -*- 



£e»^-5-Jj,m+l) 



(6.7) 



with a suitable normalization factorA/},™- This represen 



tation is useful because the states (6^) are eigenstates of 



the Hamiltonian Hq for a single chain 



Ho- — ] \j, m, k) 



|J|(l-cos(fc))|j,m,fc). (6.8) 



The computation of the matrix elements is now a bit more 
cumbersome than for the first order but follows the same 
lines of combinatorial consideratio ns. In this manner one 
finds the following counterpart of (|6.2|) for k =/= 0: 



El |j, m) = |j, m) , Ef |j, m) - |j, mil). (6.13) 

Since the matrix elements of Hn do not depend on k one 
can identify all spaces spanned by |j, to, k) for a given k 
with that spanned by |j, m). Then Hn can be written as 



Hii ^ v^ I . VU ~ m^){i + mO(j - mi){j + mi) ^ 



J' 



(hi) 



2j(2j - 1) 



V(j -m^- l)(j - mi){j + mi- l)(j + to;)' 
4j(2j - 1) 

VU + m^- l)(j + mi){j -mi- l)(j - mi) 
Ai{23 - 1) 



ErEf 



E-E^ 



(6.14) 



-ikx 



{j,m,k\S- \j,m,0) 



(j -m){j + m) 



2j 



2j-l 



(6.9) 



Still, for A2 > 1 there are excited states above the ground- 
state(s) for a given magnetization which have a gap which 
closes rapidly with increasing system size L. Therefore, de- 
pending on the chain length L, the magnetization and the nu- 
merical accuracy, groundstate degeneracies may appear to be 
threefold in the region A2 > 1. 



Here the to^ denote the quantum numbers of the state 
on which_Hr/ is acting. The combination HJjHjj appear- 
ing in ( 6.11 ) can be expressed by spin operators. In the 
isotropic case A2 = 1 a simple form can be found 



hIhii _ ^ (J, JO' - 2j(j - i)J.Ji + r'(j - 2) 



E 



4f{2j - 1)2 



(6.15) 
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which reduces for large chain length L = 2j ^ I to 
X]{i ;)(JiJ//j^ — 1)^/16 that is to the exchange between 
two chains squared instead of the simp le e xchange in first 
order perturbation theory given by (6.4). The squared 
exchange has been used as a tool to take into account 
quantum or thermal fluctuations (see ^ and references 
therein). 

So far, everything is valid for a general underlying lat- 
tice. Let us now comment on the results obtained from 
this second-order approach for three chains. For this case, 
one can easily determine the groundstate |i/'i) of ( |6.4[ ) nu- 
merically for L up to a few hundred when A2 < 1. For 
short chains, t he fi rst an d sec ond order expansion terms 
obtained from (3.4) and (6.11) compare favourably with a 



7 Exact diagonalization results for six chains 



direct diagonalization of the Hamiltonian (2.1). However, 
t he se cond-order corrections to the energy obtained from 
11) behave as L^ for large L. Alternatively, the second- 



(i 



order corrections to h{{M)) computed from these energies 
are proportional to L for long chains. This means that the 
radius of convergence shrinks rapidly to zero with increas- 
ing L. The additional factor L which limits the second- 
order weak-coupling approach for three coup led chains can 
be traced to the energy denominator ( ^.12 ) which is the 
main source of the system-size dependence. 

With hindsight this divergence of the second order 
can be expected since a one-dimensional ferromagnet can- 
not be stable to a weak antiferromagnetic perturbation if 
the chain length is long. Second order perturbation the- 
ory connects only two chains, if one compares with con- 
tributions to the ground state energy ex J one gets 
from a Holstein-Primakoff analysis for two ferromagnetic 
chains coupled antiferromagnetically one notices imme- 
diately the difficulties. Further in the limit of vanishing 
J', one would find a linear behaviour of the magnetiza- 
tion (6.5) in the thermodynmic limit. Now, if one takes 
this limit L ^ 00 &i fixed J' > 0, the asymptotic DN- 
PT square root behaviour (M) — Mc oc \/\h — hc\ which 
is characteristic for one dimension [£3|4| has to be re- 
covered. Such an abrupt change in the functional form is 
simply not compatible whith a convergent perturbation 
expansion in the thermodynmic limit. Only in t hree and 
higher dimensions, the functional behaviour (6.5) remains 
valid for the interacting quantum system (see section 7 
of [£l| and references therein). So, the second-order cor- 
rection has a chance to be convergent only if the thermo- 
dynamic limit is taken in all three directions simultane- 
ously. On the one hand, this is the case most relevant to 
CsCuCla. On the other hand, already determination of the 
first order groundstate involves the solution of a spin-L/2 
problem on the triangular lattice. Then the numerical ef- 
fort becomes much larger than in the three-chain model 
and one gains little by the perturbative treatment in com- 
parison to a direct numerical treatment of the full problem 
(which will be carried out in section 8.2). Therefore, here 
we do not pursue evaluation of the second order for the 
full three-dimensional situation further. 



Now we briefly look at the Hamiltonian (2.1) for N — 6 
chains. The coupling of the nearest neighbour pairs {i,j) 
(1 ^ hj < 6) is specified by the lines in the following 
picture: 

3 



(7.1) 




This is basically a refined version of the TV = 3 chain model 
discussed in the preceding sections and thus a natural step 
in the direction of the three-dimensional compound. In 
particular, one can see that N must be a multiple of three 
if one wishes to expect a plateau with {M) = 1/3 just on 
the basis of counting the chains. For A^ = 6, this count- 
ing argument indicates the possibility of having plateaux 
with (Af) G {0, ±1/3, ±2/3, ±1}. If one now diagonalizes 
a single layer (L = 1) of the model specified by (2T) and 
(7.1), one finds that all these values are indeed realized 
one after the other in the presence of a magnetic field. 

In a way similar to section 4 we have computed the 
complete magnetization curve of this six-chain model for 
L < 4 layers and parts of this curve for L < 35 at two 
values J'/\J\ = 1/3, 5/7 and A = 1 f\ The results for 



3 < L < 6 are shown in Fig. 7.1. The thick full lines 
in this figure are somewhat speculative extrapolations to 
L = 00 which need discussion. 

For all the sizes which we have investigated, there is 
an (M) = plateau (or equivalently a spin-gap), which 
however vanishes rapidly with increasing L. We are there- 
fore confident that there is no (M) = plateau in the 
thermodynamic limit of the six-chain model, but rather a 
steep increase arou nd (A/) = as indicated by the thick 



fuU hue in Fig. 7.1 



To see what happens to the other two plateaux which 
are permitted according to the aforementioned counting 
argument {i.e. (M) = 1/3 and (M) = 2/3) we have per- 
formed extrapolations of the available data for the bound- 
aries of these plateaux using the vanden Broeck-Schwartz 
algorithm (see e.g. ||2£|). This works better for (M) = 2/3 
where we can use data for 2 < L < 6 than for (M) = 1/3 
where just L = 2,3,4 is accessible. The result for (A/) = 
1/3 is compatible with a plateau of zero width while for 
(Af ) = 2/3 the error estimates indicate a non-vanishing 
plateau. The conclusion of the (Af ) = 1/3 plateau being 
absent in the thermodynamic limit is further supported 
by the observation that at finite size this step is about 
as wide as the neighbouring ones. On the other hand, the 



^ The average coordination in the six-chain model ( [7.l| ) is 
larger than in the three-chain model which suggests to also 
investigate smaller J' . However, the data is already somewhat 
ambiguous at J' /\J\ = 1/3. Therefore, we have added a larger 
value of J' rather than a smaller one. 
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width of the {M) = 2/3 step is notably larger than that 
of its neighbours - at least for J' /\J\ = 5/7. 



<M> 1/2 




Fig. 7.1. Magnetization curves for six coupled chains with 
Zi = 1 and a) f /\J\ = 1/3, b) J' /\J\ = 5/7. The number of 
layers is L = 3 (dotted line), L = A (thin full line), L = 5 
(dashed line) and L = & (dashed-dotted line). The thick full 
line is the extrapolation to L = co discussed in the text. 



Therefore we have decided to draw plateaux at (M) = 
2/3 in Fig. 7.1 but none for (M) = 1/3. The boundaries 
of the former case are the only points where fairly reli- 
able extrapolations to i = oo are available. Otherwise we 
follow the procedure of fSy,^ and draw a line through 
the midpoints of the steps at the large st a vailable system 
size, although it is apparent from Fig. 7.1 that finite-size 
effects are still of some importance. 



Finally, the transition to saturation is of the same type 
as that for the three-leg ladder, i.e. it is again compatible 
with the DN-PT universality class |g^,|j, but subject to 
strong crossover effects in the region J < 0. 

The main result of this discussion of A^ = 6 chains 
is that already doubling the number of chains leads to a 
vanishing of the {M) = 1/3 plateau - at least for Z\ = 1 
and weak coupling J'. This is a first indication that this 
plateau which is easily observable in A^ = 3 coupled chains 
might actually not survive the infinite-lattice limit in the 
plane perpendicular to the chains. 



8 The three-dimensional case 

The results from the various techniques which we have 
applied so far to simple models agree well with each other. 
This encourages us to proceed a bit further and look at 
the three-dimensional model correspondi ng t o CsCuCla. 
This means we consider the Hamiltonian (2.1) where the 
nearest neighbour pairs (i, j) will now be those of the two- 
dimensional triangular lattice. 



8.1 Ising expansions 

We start with Ising expansions. The structure of the state 
with {M) — 1/3 and the lowest single-spin excitations 
above it is readily inferred from the treatment of three 
coupled chains and the two-dimensional triangular lattice 
[ pl[ . It is then straightforward to obtain fourth-order series 
for these excitation energies^. The explicit series can be 
found in appendix A. As before, these results correspond 
to the boundaries of the (M) ~ 1/3 plateau for sufficiently 
small A. At large values of Z\, the associated transitions 
should be first order such that one would have to look at 
other excitations to determine these boundaries. 

If we inse rt the va lue J' = \J\/?> into the raw fourth- 
order series (A_^,A^) we find that the plateau closes at 
Ac — 0.939. Using instead J' = | J|/6 (which is appropri- 
ate for CsCuCla), the ending point shifts to Ac = 0.948. 
Both values are slightly smaller than those obtained for 
the three-chain model with the same order. However, in 
this region the actual values should not be taken too seri- 
ously, i.e. higher orders should still important and numer- 
ical data will indicate that the true Ac is probably larger. 
The main conclusion obtained from the Ising series is that 
the (M) = 1/3 plateau closes in the region Z\c w 1, proba- 
bly somewhere slightly below one in the region of XY-type 
anisotropies. 



8.2 Exact diagonalization 

To locate the ending point more accurately and to see 
if a jump in the magnetization curve develops, we have 
also numerically computed magnetization curves on three- 
dimensional 3x3 x_L clusters. Some magnetization curves 
are shown in Fig. ^.l| . Two criteria can be used to discuss 
the magnetization process in the region (Af ) = 1/3; 

1. If the magnetization curve passes smoothly through 
(M) = 1/3, the width of the finite-size plateau at 
(A/) = 1/3 should scale as 1/L (with L the number 
of stacked segments of the triangular lattice). 

2. For a smooth curve at magnetization (Af), the neigh- 
bouring steps of the finite-size magnetization curves 
should have the same width. If a plateau is present, 
the step with magnetization (M) should be broader 
than its neighbours. A jump at sufficiently large sys- 
tem sizes should in contrast be indicated by a step that 
is smaller than its neighbours at small system sizes. 

A remark is in order before we apply these criteria to 
our data. Although we expect to obtain a reasonable ap- 
proximation to the thermodynamic limit despite the small 
linear size of the systems accessible to us, one can clearly 
not rule out that the behaviour would change on larger 
systems. Our conclusions should therefore be taken as a 
probable possibility. 

For J' /\J\ — 1/6 and A = 1, both criteria indicate the 
presence of a small plateau with (A/) = 1/3. Instead for 

® This can probably be pushed a bit further. Higher-order 
series have been computed numerically e.g. for the Heisenberg 
model on three-dimensional cubic lattices fB2f. 
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order spinwave theory. The finite-size m agn etization curves 

Ic). Again, both 



Fig. 8.1. Magnetization curves of 3 x 3 x L clusters with a) 
J'/\J\ = 1/6, Zi = 1; b) J'/\J\ = 1/6, A = 0.97; c) J' /\J\ = 
0.172882193, Z\i = 0.98789825, A2 = 1. The number of layers 
is L = 2 (dashed line), L = 3 (thin full line) and L = 4 (dotted 
line). The bold full line is an extrapolation to L = 00. 



J' /\J\ = 1/6 and A = 0.97 both criteria apphed to L = 2 
and L = 3 point to a smooth magnetization curve in the 
region (M) = 1/3. This suggests I > Ac > 0.97 which is 
slightly larger than that found from the fourth-order Ising 
expansions. 

Since we did not find evidence for a jump so far, we also 
performed diagonalizations for the parameters J'/\J\ — 
0.172882193, Z\i = 0.98789825, A2 = 1, since for these 
precise parameters a jump was found in ||^ using first- 



for these parameters are shown in Fig. 
criteria discussed above are compatible with a smooth 
curve at (M) = 1/3, i.e. neither a plateau nor a jump 
seems to be present. However, on a system with volume 
27 = 3 X 3 X 3, we would expect some indication for a 
jump of size S{M) w 0.04 as estimated in Q or even 
5{M) « 0.01 as observed experimentally Q. One may 
attribute the absence of such signals for a jump to finite- 
size effects beyond the discretization of the magnetization 
curve at a volume with 27 spins. However, it can be seen 
at larger (M) that increasing the number of layers L leads 
only to small finite-size effects. On the other hand, also the 
approach of g is based on several approximations and is 
therefore not guaranteed to describe the situation with 
spin 1/2 accurately. 

To conclude this section, let us briefly comment on 
the utility of order parameters which one might suspect 
to be useful for characterizing the transition. First, one 
could check nunrerically that the correct spin structures 
have been used in g]. However, we do believe that the 
spin structures are indeed the umbrella-type and copla- 
nar configurations in the low- and high-field regions, re- 
spectively. In fact, this has been reliably established via 
a mapping to a gas of hard-core bosons 1 33 for magnetic 
fields close to the saturation value. The question which we 
raised above regards the nature of the transition which oc- 
curs at intermediate fields, i.e. if it is second order, first 
order or occurs e.g. via an intermediate phase. However, 
we do not expect much additional insight into the nature 
of the transition from a numerical determination of or- 
der parameters, mainly because of finite-size effects such 
as the discretization of the magnetization axis. For this 
purpose the magnetization (M) (which we have discussed 
above) might even be somewhat better suited although it 
is not an order parameter. 



9 Conclusions 

We have theoretically observed a plateau with (M) = 
1/3 in a frustrated triangular magnet. For applications 
to CsCuCls, it is important to observe that the width of 
this plateau depends on several factors. 

Consider for example a single triangle. If there are only 
three spin-1/2 spins, the magnetization can obviously only 
have the values (Af) = ±1/3 and ±1. The plateau with 
(M) = 1/3 survives (at least for A = 1) the ferromag- 
netic stacking of infinitely many triangles, but the inter- 
action reduces its width. A similar phenomenon occurs if 
one instead takes the limit of a two-dimensional triangu- 
lar lattice in the plane. Also here, the plateau survives the 
thermodynamic limit for Z\ w 1 |pj,El|, but again gets 
narrower. 

A further factor is the XXZ-anisotropy A. For a mag- 
netic field along the c-axis, the plateau with (M) = 1/3 
survives only a small amount of XY-like anisotropy be- 
fore it disappears. The same observation has alre ady been 
made for the triangular lattice antiferromagnet |^,^ . If 
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the magnetic field is applied in the triangular plane per- 
pendicular to the c-axis, an opposite behaviour is known 
from the triangular lattice antifcrromagnet |35 : Here one 
observes a clear plateau-like feature in the region {M) = 
1/3 even in the XY case A = 0. 

If one combines these observations with the parameters 
which are appropriate for the description of CsCuCls , one 
understands the absence of a plateau for a magnetic field 
in the direction of the c-axis and that a plateau- like feature 
only occurs if the magnetic field is applied in the plane g . 

In contrast, we failed to reproduce the jump observed 
experimentally |2| for a magnetic field along the c-axis 
with the model (2.1) in the appropriate parameter region. 
While we observe several first-order transitions for Z\ > 1, 
we have no evidence for first-order transitions for A < I. 
Actually, our results for the three-chain model as well as 
similar ones for the two-dimensional triangular lattice |gl[ 
provide quite strong evidence against such a jump for spin 
1/2 and A < 1. However, a first-order spinwave analysis 
P] predicted a ju mp in the magnetization curve of the 
Hamiltonian ( pjl ) with zii < 1. Since it may be possible 
that this jump cannot occur in less than three dimensions, 
we have also investigated the three-dimensional spin-1/2 
model numerically. Using precisely the same parameters as 
in 0, we did not find any evidence for a jump on 3 x 3 x L 
clusters either. This discrepancy is puzzling, in particular 
since the weak-coupling analysis shows that weakly cou- 
pled ferromagnetic chains are close to a classical situation 
(though in two rather than in three dimensionsll) . While 
we cannot fully exclude that one would find a jump in the 
spin-1/2 model at larger system sizes, also the discussion 
in Q is based on several approximations. We therefore be- 
lieve that the issue of the experimentally observed jump 
deserves further theoretical attention. One possibility is 
also that this jump is cause d by effects which are not in- 
corporated into the model ( ^.l|) , such as elastic response 
to the external magnetic field. The anomalies observed in 
a recent soundwave propagation experiment [p6[ may be 
a further indication for the necessity of such extensions. 

In summary, we have shown that a three-chain model 
describes the low-temperature magnetization process of 
CsCuCla almost as well as the full three-dimensional model. 
The three-chain model has the advantage of being sim- 
pler to analyze. We have concentrated on zero (or small) 
temperature, but there is also a large amount of experi- 
mental data on the h-T phase diagram (see p7| , p8||3S( | for 
a small selection of recent results) which are waiting to 
be discussed in the framework of microscopic models. The 
three-chain model which we have discussed in the present 
paper may be a useful starting point since one should be 
able to treat it with a finite-temperature variant p0| , ^ , p2[ 
of the DMRG procedure. 

From a technical point of view, we have shown that 
the Hamiltonian (2.1) with spin 1/2 can be analyzed by 
several methods. Each of them has the potential of be- 



'^ More precisely, each ferromagnetic chain can be regarded 
as one single spin with large S. It might therefore actually be 
interesting to try to treat our second-order expression in J' in 
three dimensions using spinwave techniques. 



ing pushed further to obtain more accurate results. This 
applies in particular to the weak-coupling approach: The 
second order diverges for the three-chain model as L — > oo, 
but it remains to be investigated if it converges in three 
dimensions. The computation of other quantities than just 
the magnetization is also straightforward with each of the 
methods used in the present paper. 

Useful discussions with D.C. Cabra, T. Nikuni, I. Peschel, 
P. Pujol, U. Schotte, H. Tanaka and M.E. Zhitomirsky are 
gratefully acknowledged. We are indebted to the Max-Planck- 
Institut fiir Mathematik, Bonn-Beuel and the C4 cluster of the 
ETH for allocation of CPU time. 



A Ising series 

In this appendix we present explicit series in the Ising 
anisotropy A~^. 

For the (Af) = 1/3 plateau in the three- chain model, 
one finds the following sixth-order series for the gap of a 
single spin flipped up 

'J 1 ^-1 , 5J'-2J^_2 



h,,^A{J'-J) + J + J'U—A 



16./2 



164 JV^ -~ 258 JV' + 2767-^ - 93J'^J + 19.1''^ 3 
192 J3 (J' - J) (J' - 2 J) 

/\-4 , 

5 7 (19050 J' J"* 

2304 J4 (J' - jy (J' - 2jy V 

-82347'" J3 + 3297.1''" J^ - 27548J'V^ 

+20232 J' J*^ - 1136 J'V + 171J''^ - 6480 J^ 

A-^ 

55296 (J' - 2jf J5 (J' - 3 J) (J' - 4J) (J' - jf 

(62217652 J'V^ - 35026088 J'V^ 

-76513096J'V^ + 65015952 J'V^'' 
+ 14995748J'V^ - 7562886J'V^ 
+8581J'^^ 4- 4998735 J'V" - 25546107'^^ 
+7833807'^" J2 - 128008J'" J 

-35276256 J" J' + 9134208 Ji^) i 

+OiA-^). 



(A.l) 



For an excitation with a single spin flipped down one finds 
instead the following sixth-order series 

^ T T J T/2 2 J — 3 J . 9 



+ 



47.PJ'^ - 39JV -I- 30J'^ - 25J'^ J + 5J'^ 
64J3 (J' - J) (J' - 2 J) 

A —A 

/3 ^^ 



A' 



4608 J4 (J' - jf (J' - 2jf 



24218 J' J 
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-177767'^'"^ + lOgOlJ'V^ _ 30028J'V^ 
+24264J'J^ - 4054J'V + 611J'^ - 8784J^ 

Ar^ 

110592 (J' - jf (J' - 4J) (J' - 3J) J5 



1 



(J' - 2 J)^ (2 J' - 3 J) 



( 468904952 J'^ J® 



-400904412J'V^ 



488731468J'V^ 



-415953144J'V^° - 263997936J'V" 



-309762362J'V^ 



199502115J'V^ 



-23421312J 



13 



51214J' 



,13 



111406752J1V 



+97003020 J'° J"* - 32952968J'^" J3 
+7259020J'" J2 - 923565 J'^V") i 

+0{A-^). (A.2) 

In a similar way, one obtains the following fourth-order 
series for the single-spin excitation above (M) = 1/3 in the 
three-dimensional model 



J A 



J 



3J^ 
2 
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8A (J' - J) (3 J' - 2 J) (2 J' - J) (J' - 2 J) 

/2 
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-912414730048J'V" + 1875316243840J'V^" 
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(A.4) 



(A.3) 



and 



/lc2 = (3 J' - J) Z\ + J - 



3J'^ (6/2 - IIJJ' -I- lOJ'^ 
AA [J' - J) (2 J' - J) (J' - 2 J) 



The results ( |A.3 , |A.4 ) are consistent with earlier ones. Firstly, 
for J = 0, one recovers the series for the tri angu lar lat- 
tice ML For ge nera l J < also (A.l) and (A.2) match 
with ( [A.4|) and ( [A.3| ), respectively up to first order in J' 
(one just has to rescale J' by a factor 3 to account for the 
different coordination number). 



B Details of the DMRG procedure 

Following the standard infinite-size algorithm of the DMRG 
method |0| , we enlarge the chain in each iteration step by 
two triangles (in total six spin-1/2 sites), i.e. by 8 x 8 = 64 
states. This is different from ^^ wh ere t he totally antifer- 
romagnetic case of the Hamiltonian {2X) has been treated. 
There, the chain was built up by adding a total of only two 
spin-1/2 sites in each step, so that it took three DMRG 
steps to add two complete triangles to the original lad- 
der. In this procedure more truncation operations of the 
Hilbert space had to be performed than in our calculation. 
But by adding eight-state sites in each step the number of 
basis states discarded in the truncation process is larger 
than in the case where only two-state sites are added. As 
far as we know it has not yet been tested which procedure 
is more favourable. 

In order to limit the computational effort, we applied 
only the infinite-system algorithm even though more ac- 
curate results could probably be obtained when one would 
use in addition the finite-system algorithm. To calculate 
the boundaries of the magnetization plateau we have used 
spin conservation and performed three different calcula- 
tions to target each of the groundstates with magnetiza- 
tion {M) = 1/3-1/S'tot, 1/3, 1/3-^1/5*401 separately. The 
difference between the groundstate energy in the sector 
{M) — 1/3 and the energies in the neighbouring sectors 
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then gives the magnetic field: 
/ici = Ea ((M) = i + ^, i) - i^o ((Af ) = i, L) (B.l) 
hc2 - £^0 {{M) = \,L)-E„ Um) = 1 - ^,l) .(B.2) 



Fig. p3.1| shows a typical spectrum of the eigenvalues 
of the density matrix for such a calculation. 
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Fig. B.l. Eigenvalue spectrum of the reduced density matrix 
of a ladder with 24 spins (L = 8), calculated with 64 kept 
states in the subspace (M) = 1/3. The couplings are J = —1, 
J' = 1/3 and Z\ = 1 



The eigenvalues drop exponentially as is known from 
other quantum spin chainsn and thus gives rise to a very 
small truncation error. On the other hand, we found that 
at least 160 states of the density matrix, i.e. 1280 states 
in one part of the chain had to be kept to obtain reliable 
results up to three decimal places in the energy. The eigen- 
values of the density matrix lie about three to four orders 
of magnitude above the eigenvalues in the corresponding 
spectrum [|6| of the totally antiferromagnetic case studied 
in |26|| . This is one reason for the lower accuracy of our 
DMRG calculations. There is also a big difference between 
the truncation error and the real error in the groundstate 
energy which indicates that the structure of the ground- 
state could be approximated better by applying the finite- 
system algorithm. 
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